“Calhoun 


Institutional Archive of the Naval Postgraduate School 





Calhoun: The NPS Institutional Archive 
DSpace Repository 


Theses and Dissertations 1. Thesis and Dissertation Collection, all items 


1967 


A method for evaluating seakeeping 
performance of destroyer feasibility designs. 


Kalb, David G. 


Webb Institute of Naval Architecture 


http://ndl.handle.net/10945/12653 


Downloaded from NPS Archive: Calhoun 


: Calhoun is the Naval Postgraduate School's public access digital repository for 
/ (8 D U DLEY research materials and institutional publications created by the NPS community. 
«ist : Calhoun is named for Professor of Mathematics Guy K. Calhoun, NPS's first 





KNOX appointed — and published — scholarly author. 


i LIBRARY Dudley Knox Library / Naval Postgraduate School 


411 Dyer Road / 1 University Circle 
Monterey, California USA 93943 





http://www.nps.edu/library 


aed tal 
one “tite 
ens 
t payt h 


ah 


NPS ARCHIVE 


; i 1967 
Seed 8 =6KALB, D. 





A. METHOD FOR EVALUATING 


SEAKEEPING PERFORMANCE OF | 
DESTROYER FEASLBLLiTY DESIGNS 





LT... DAVID G. KALB 1467 


Thesis 
K115 








DUDLEY KNOX UBRARY 
NAVAL POSTGRADUATE SCHOOL 
MONTE Rey ca, 93943-5104 








_ A METHOD FOR EVALUATING SEAKEEPING 






“PERFORMANCE OF DESTROYER PEASIBILITY DESIGNS 
A THESIS 
SUBMITTED TO THE FACULTY OF 
WEBB INSTITUTE OF NAVAL ARCHITECTURE 
IN PARTIAL FULFILLMENT OF THE REQUIREMENTS FOR THE 
DEGREE OF 
MASTER OF SCIENCE 
IN 
NAVAL ARCHITECTURE 
by 
Le DAVID G. KALB, USN 


30 MAY 1967 





WES A€cyiVve 
ie 


ALE, O- 





LIBRARY 


NAVAL POSTGRADUATE SCH | 
OOL 
MONTEREY; CALIF. 93940 | 


ACKNOWLEDGEMENTS 


The author would like to express his sincere 


appreciation to: 


Captain Robert A. Hinners, USN (Ret:), who, as thesis 
advisor, Was always available for advice and encouragement; 











Professor Edward V. Lewis, thesis advisor, for providing 
technical advice and a generous amount of his time; 


Mr. Charles Maiorana, DTMB, without whose help the pro- 
gramming, trouble-shooting and numerous computer runs could 
not have been carried out; 


Members of the Research Department of Webb Institute for 
their advice and assistance in analysing the computer 
printout. 


Manaliy, to Alice. 





TABLE OF CONTENTS 


Abstract 
Nomenclature 

List. of Figurés 
tntrodwerion 

The Problem 
Solution-Part I 
Solution-Part II 
Solution-Part III 
Analysis of Results 
Conclusion 
Discussion And Recommendations 


References and Bibliograpny. 


Appendices ° 


A Sea Spectrum 
B Input for Seakeeping Program 


C Sample Calculation Demonstration 


Confidential Annex - Separate Cover 


10 
14 
ty, 
Be 
23 
2/7 


42 
LU 
4? 





ABSTRACT 


The author discusses a method of determining ship 
motion responses for a number of destroyer-type hulls which 
can be evaluated, relative to each other, for their per- 
formance in a given sea state. The method is intended to 
be incorporated in and be compatible with a computer-aided 
destroyer feasibility study. 

A program of the method, along with a discussion, 
and sample results are presented. Output from the program 
using the Neumann and the [.T.1T.C. Wave spectra are con- 
pared. Results and supporting data are presented and 


discussed. 
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INTRODUCTION 


In 1964, when long-range studies of destroyer systems 
were begun at the Center for Naval Analysis (CNA), it was 
realized that many systems and their capabilities would have 
to be considered and analysed. In order to support lonz- 
range analyses of costs, effectiveness, trade-offs and sen- 
sitiveness of alternative designs, principal characteristics 
of these ships Were required. There were two methods of pro- 
viding these data: manually, which limits the amount of 
data on a cost or time limitation, and by computer, if a 
mathematical model could be developed suitable for computer 
computation. 

Since a large number of ships was to be investigated, 
the manual method was considered impractical. CNA, with 
the full cooperation of the Conceptual Design Division (CDD) 
of the Naval Ship Systems Command (NAVSHIPS), undertook to 
develop a suitable parametric design model which closely 
followed the current U.S.Navy design procedures ’and con- 
struction practices. From this work evolved the present 
"Destroyer Feasibility Study" program currently in use at 
NAVSHIPS. 

The problem of developing a computer feasibility 
study is one of determining a set of principal dimensions 
and form coefficients which simultaneously satisfy oper- 


ational, propulsion, Weapons, electronics, weight, moment 





2 
and volume requirements. The solution to the problem in- 
volves the consideration by computer of a wide range of 
combinations of these ship Ee ee A unique solu- 
tion does not exist for a given set of military requirements, 
rather, a family of ships can be determined for a spectrum 
of hull lengths and form coefficients. 

The computer feasibility study program is a design 
tool which allows the analyst-e-designer to consider a varied 
and wide assortment of ships which are all compatible with 
the initial conditions of the input data. The designer 
then must take these ships and narrow the spectrum for a 
more detailed study. He arrives at his selection and elim- 
ination on the basis of costs, powering, payload (weapons) 
and naval architectural considerations based on current 
requirements, cost and effectiveness studies, mission ful- 
fillment and his knowledge and experience. He then has 
a limited group of ships from which the final design may 
evolve. 

The present NAVSHIPS program does not utilize any 
seakeeping performance data. It would help the analyst-designer 
if seakeeping performance information were supplied along with 
the other output data so that in his determination he could 
€évaluate the operational potential of the various ships under 
consideration when compared to each other in the various sea 


states. 











THE PROBLEM 


The output from the aforementioned study lists a number 
of sets of design parameters from which a series of ships 
could be designed, each compatible with the requirements 
used as inputs. (Table I) However, no lines are drawn, 
there are no body plans, and hence the ships have no delin-~ 
eated form. Seakeeping programs currently in use utilize 
the strip theory of integration which requires that the ships 
undergoing analysis have a body plan or at least a water line, 
half breadths, draft and areas at stations. 

The first part of the problem, then, is to develop 
a program which will give the feasibility ships either a hull 
form or provide sufficient data from the dimensions and co- 
efficients which then can be utilized in a seakeseping study. 
The second part is to utilize the hull form developed as an 
input to a seakeeping analysis computer program. The third 
part is to-analyse the results of the seakeeping output 
and define useful indices or criteria which will help the 


analyst-designer in making his final choice, 
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A solution to the first part of the problem involves 
the development of a method which will compute the hull 
form inputs (beam, draft and sees area coefficient for 
each of 20 stations along the hull) required for a sea- 

Keeping program. 

A current project at NAVSHIPS is the updating and 
revision of the feasibility study. program. From a comprehensive 
study of the destroyer-type ship conducted over a period 
of years, a parent hull was developed by Mr. James Milis!, 
From this. parent, four other hulls were derived by Mills, which took 
into account the upper and lower ranges of destroyer-type 
parameters. The five ships were then placed in the Y-G 
Hull program 2 and relationships among the hulls were 
analysed for the purpose of refining deck area and volume 
Calculations in the feasibility study. 

This previous utilization by Mills of a parent and four 
derivatives for the non-seakeeping part of the feasibility 


study suggested a possible modus operandi to the author for 


1 see Confidential Annex. 


The Y-G Hull computer-program was developed by G. Hirschberg 
and Lt. R. Warters, for computing any or all of six options: 
1) Hydrostatics and Bonjean curves, 2) Compartment vol- 
umetric properties, 3) Longitudinal Strensth, 4) Floodable 
Memetn, 5) Limitins Drafts, 6) Intact cross curves of $Sta- 
Drlity. The primary input data are the offsets of the hull 
and LBP. 








arriving at a method of computing the needed hull form or 
input parameters when lines are not available. An empirical 
relationship from this series for area, beam and draft seemed 
possible and this was the approach taken. 

The two parameters which had been varied in the five 
ships were prismatic coefficient and midship section coeffi- 
cient. In order to avoid a three-dimensional problem, one 
of the variables had to be held constant to reduce the prob- 
lem, to two dimensions. 

After a series of discussions with personnel in the 
Hydromechanics Laboratory of David Taylor Model Basin (DTMB) 
concerning the effects of Cp and Cx on ship motions, 
it was decided that the effects of varying Cx were very 
small when compared to varying Cp.: Cx was then chosen 
as the one to be kept constant for the purpose of develop- 
ment of the interpolation program. However, it was anti- 
Cipated that a variable Cx could be reincorporated in 
the final program. 

With due consideration for the above decision the 
following method was developed. 

Values of local beam and draft for the twenty stations 
of all five ships were lifted from the body plans and non- 
dimensionalized by dividing by Bx and Hx, Areas at each 
station were computed by the Y-G Hull program and also non-~ 
dimensionalized by dividing by Ax. These non-dimensional 
values were then plotted against Cp, the variable, for each 


of the 20 stations (Figs. C-1, c-2, @e3).2 


3Confidential Annex 





Simultaneously ttheseivalues were put into.a 

least squares polynomial curve-fitting computer program 

and equations using Cp as the variable were computed for 
each of the twenty stations. Twenty equations for non- 
dimensional area, twenty equations for non-dimensional beam 
and seven equations for non-dimensional draft were computed. 
(Draft is constant from station zero to station 14 where 

the transom line intersects the base line.) 

The resulting equations were checked for accuracy 
and then written as a controlling program for supplying 
hull form inputs to a seakeeping program (to be discussed 
in the next section) currently available at DIMB. 

The controlling program solves the series of equations 
using Cp as the variable. The resulting non-dimensional 
values of area, beam and draft are then dimensionalized 
by multiplying by the actual input values of maximum section-~ 
al area, maximum beam and maximum draft. Section area co- 
efficients are then computed for the twenty stations. The 
computed values of beam, draft and section area coefficient 
are then stored in the computer-memory for use by the sea- 
keeping progran. 

At this point the question of the validity of rein- 
troducing Cx as avariable in the equation: Ax = Cx . Bx. Hx , 
arose. (It will be recalled that Cx was temporarily kept 
constant for the purpose of developing the interpolation 


curves and resulting equations). To justify the validity 





of this reinireduction, chee pananéitere of two of the five 
original ships, which varied from the assumed constant 
value of Cx, were run through the control program. A 
comparison of this output with respect to sectional areas 
Was made with the more exact areas computed by the more 
detailed Y-G Hull calculations. The differences were very 
small (See confidential annex) and are insignificant when 
considering ship motions. This comparison led to the con- 
clusion that Cx can be reintroduced into the program as 

a variable when computing maximum section area. 

The first part of the original problem was then 
considered solved. The input variables of Cp, Cx, Bx 
and Hx could now be developed into a mathematical hull which 
Was compatible with the ships arrived at in the feasibility 
study, using the output of the feasibility program as input. 
Section area curves as computed by the control program are 


shown in Figures [I and II. 





THE SOLUTION — PART II 


The seakeeping program selected for use with the 
control program was developed at Massachusetts Institute 
of Technology under a NAVSHIPS contract as part of the 
Fundamental Hydromechanics Research Program administered 
by DTMB.(18] It was originally written by K. Haslum, making 
use of subroutines from Stevens Institute of Technology 
Which utilize the Korvin-Kroukovsky theory for computing 
added mass and hull damping coeffigients. It was later 
expanded by L. Vassilopoulos to include statistical re- 
sponses in irregular seas applying a Neumann energy spectrun. 
This program will be referred to hereafter as the MIT pro- 
sram. This program was revised by N. Salvesen (11) at 
DIMB, to include a plot routine (not utilized in this study), 
changes in input-output format, changing to FORTRAN IV com- 
puter language and adding speed contributions to the com- 
puted coefficients according to the theory of Tuck. (The 
program still remains basically the MIT program). 

The MIT Program computes pitch and heave amplitudes 
and their phase angles in regular head seas, calculates 
bending moments and computes several statistical responses 
in irregular head seas. The particular statistical response 


used in the present study was the root mean square value 


AUnpublished work by Tuck at DTMB 
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of relative vertical motion at station 2. The relative 
wertical motion réfers to the distances bety- en Gena er 
surface and the DWL as the ship responds to the seaway. 

The MIT Program was used for determining the statis- 
tical responses for thirteen (13) hull forms; the original 
five ships from which the equations previously mentioned 
were derived plus an additional eight (8) test ships (TABLE 
II) taken from the output of the Destroyer Feasibility 
program, mentioned in the Introduction. 

The author then changed the subroutine that calcu- 
lates the sea spectrum by the Neumann equation to a new 
subroutine that utilizes the International Towing Tank 
Conference, 1966 [22] , equation, changing from wind speed 
input to wave height input. (FigurelII) The 13 ships were 
then run through this program. 

The input information to both programs mentioned 
above is the same with the exception that significant Wave 
height values are substituted in the latter for significant 
Mind speed values in the former. The ITTC spectra 
were used by the author since they are more realistic. A 
comparison of the two spectra and the resulting response 
spectr@® is shown in Figure IV,for one case. 

Information not calculated by the feasibility progran, 
radius of gyration and longitudinal center of gravity, was 


set at .25L and .50L respectively. These values are commonly 
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used in initial design studies. Other input values are 


as follows: 


Number of stations 20 

Ship speed _ 10 to 30 knots 
Significant Wave Height (Hy ) 14.5 to 48.0 feet 
Frequency of encounter * 10 to 3.10 rad./sec = 


These were held constant for all ships. The variables were 
the ships themselves with beam, draft and station sectional 
area coefficients computed by the control program described 


in Solution Part I. 


*Calculated range to encompass the maximum area under the re- 
sponse spectrum which would give accurate values of response. 





Bh SOLUTION - PARP as 


The format of the printout of the seakeeping pro- 
gram, described in Appendix A, allows the analyst to select 
from a number of statistical values of response. The square 
root of the area under the squared amplitude density spectrum 
of relative motion at the "slamming station" (.1L aft of FP 
or station 2) is calculated by the program. These r.ms. 
values were tabulated for each ship and used in line 3 of 
the probability tabulation. (TABLE V) 

The seakeeping criterion chosen for investigation 
is sonar dome immersion for a dome mounted at .1L from 
FP at the ship's base line. 

If one accepts the hypothesis that the sonar will 
always function properly if its dome remains at least one 
foot below the water surface, a criterion of performance 
can be established which expresses the limit of operability 
of bow mounted sonar in terms of the frequency that the 
relative bow motion exceeds the draft of the vessel minus ; 
Ome root. If this criterion is applied directly to the cal- 
culated r.m.se motions, a limiting r.m.s. value can be selected 
which has a value of (H - 1) feet for each ship considered. 

A Wave height frequency of distribution of the spectra 
can then be selected. The distribution used by the author 
is based on that given in [2] for the North Atlantic and 


is given in TABLE IV. This is not a ‘particularly severe 
{ 


14 








7 fF 

















4dd°= B OO 0°0 Lcz* Z 
020 0°0 00°T 9 
Sqouy o€ S0L°=9729°02- —=Ss— OL * LZ TI- 60° S 
paeds Loh = 2 ee 9£0° Tae 1 
dyus 0S°02 GT°QT OO°E € 
¢00° 0z00° ded° a 
608°= 0°O OmO LGz° d 
0°O 0°0 olen § 9 
sjouy Sz nS9°=b9n°0Z= 65°9T= ne? 4 
pseds 680° t= Oo: t£O° 90°2 ty 
dtus 8z°61 C2°9T og°2 C 
¢00° 0z00° Geue z 
O°Sh ET Th ai : 
OSe*=28 0°0 0°0 Led° ra 
070 0°) evel | 9 
Suet O¢ 009°=2 S8"0¢c- GIO 1946 G 
' paeds 6nd°c=Z 600° 150% 4g°t qT 
dtus T° Qt 65°ST tome Con 
G0O0° 0z00° A Sf G 
0° 8t Cay 7 i Gent T 
ee = A Men 9 eee OK UUULULU!.ULULU.!DLUL.ULULULLULULUUE 
gZg°=s 0°0 0°0 : : ° Lous d 
0°O 0°O 00° 9 
eC Ue SI 60S °=© S9°1e= GCs g9°d S 
peeds Cer°€=S 600° Of 0 89°T t 
dyus 06 29 °HT Soe € 
S00° Waleroy STi 2 
O°St gy C°Ht T 
069°= 2 Beea 0 Led° 9xzseutT 2 
Ome 0°O 00°T 6°CG=SHI SUpPpesoxe 
qou jo AATTIQeqorIg g 
21S°=) L6°¢oae Gee Ges 7 oa a) 
ren 66T°€= << 600° 620° Bout €xz SeUuyt yonporg + 
Lt ae 99°2T L547? 7 Ql°c UOTIOW* TOI enTBaA sBI € 
ae S00° OiGD.0 * SG Fo Aouenbety eAumM Z 
0° St C° id C*Ht *VJ-£I0G94BD OABM YT 


(47 out a nati 


XN TA TTIW TO. GT 


























iy 


distribution and could easily be altered as desired to 
represent other ocean areas or heavy weather conditions. 

To perform a probability analysis of the ship's 
motions a measure of the variability of r.m.s. response 
mithin each sea condition 1s required. A value of the 
ratio of standard deviation (ow) of r.m.s. relative motion 
to mean r.m.s. value equal to 0.16 was selected from [2]. 
To demonstrate the method and reduce calculations this 
value of S/m was used for all ships at all speeds. More 
accurate V/m values can be calculated DY COMputing Ship 
responses to a randomly selected group of sea spectra in 
several ranges of sea severity as was done in[{2]. 

It Was assumed that the destroyers being analysed 
could perform their mission if their sonar Was operable 
for 964 of the motion cycles over a particular short term 
Sieuetion, for example, one watch during a patrol. This 
means that 4% of the relative bow motion peaks oy her tae 
would be greater than (H - 1) feet, causing the sonar dome 
to rise to less than one foot below the surface, making it 
inoperable. 

If relative bow motion is Rayleizh distributed and 
subject to a limit of (H - 1) feet with a 96% probability 
of Operable sonar time, the corresponding limiting r.m.s. 


value from probability theory is ,39 times (H - 1). In the 


case shown in TABLE V, where draft is 16.31 feet, 15.31 x 


0.39 = 5a. 
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Taking all of the above into account, a tabular for- 
mat for calculating the probable performance of the snips 
was set up as given in [23] and used in (2) » Five ships 
with the widest variance of parameters were used, in lieu 
of calculating all eisht, in order to show clearly the dif- 
ferences in performance. The result of each tabular calculation 
Mme the overall percent of time the ship could be expected to 
perform its mission in different head seas in the assumed long 
term wave distribution. <A sample tabulation is shown in TABLE 
V for one particular ship at several speeds. 

The results of all the tabulations are then plotted 
as a function of ship speed (Figures V & VI) and show the 
overall performance of each ship over the speed range. The 
figure shows that desizns 5 and 8 have definite superiority 
and could patrol in the North Atlantic at speeds above 25 
knots better than 77% of the time while achieving satisfactory 
sonar performance. 

The method proposed by the author results in a pre- 
diction of the percent of time a ship could be expected 
to perform a given mission satisfying a particular seakeeping 
Sereerion. <A group of ships with varying characteristics, 
such as the feasibility prozram provides, can thus be compared 
relative to each otner based upon their predicted performance. 

The criterion used in this study is an example, and 
other limiting values and probabilities may be more anppro- 
priate. The r.m.s. value used in line 5 of TABLE V can be read 


from Figure VII for any limiting value and any probability. 





ANALYSIS OF RESULTS 


The final results of Chis stu@y are tne pert ormance 
curves which show the predicted performance of the shins. 
The trends of performance for the variation of hull char- 
acteristics are consistent with trends studied in [19}, [23] 
and [2] » The longest ships, 5 and 8, are shown to per- 
form the best of the eitht feasibility ships. Since dimen- 
sions and coefficients, other than length, are very close 
in these desicns, tnis result is believed to represent 
the effect of length. Ship number 5 is slirhtly better than 
8, on the performance curves. (It will be noted that on 
FIGURES VIII thru XIII ship number 5 shows a greater relative 
motion than snip number 8. This would lead one to believe 
that ship number 8 is the better ship. However, for the criterion 
Meger invesctization ship number 5 performs better. This is due 
bO the deeper draft of number 5. The deeper draft then, offsets 
j@@emclizht disadvantage of motion for this criterion.). Desien 
Number 1 is the worst performer due to a number of possible 
fj tevors: amonz them, shorter length, staller Cp and reduced 
beam to draft ratio. 

The author feels that the variation in length amone 
the 8 feasibility designs overrides the effects of small 
Variations in otner parameters. Here is where the analyst 
can apply nis knowledze. Tne lLonzer ships may perform 
better but they usually cost more. Is the increased effect- 
iveness of designs 5 and 8 worth the cost of increasing 
lenzth? Using Figure V, a 4% increase in effectiveness 
of these designs over desizsn 1, at 25 knots may not yYustify 
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the increase in length of 25 feet. That decision is left 
to the analyst. 

For the five JLM ships, length, beam and draft 
were constant thereby allowing a check on proper trends 
with the variation of Cp and Cx alone. Trends Nave 
been studied by L. Vassilopoulos [19] for a series of 
Mariner-type hulls. This study showed an increase in Cp, 
While reducing Cx reduces ship response, and an increase 
in Cx while reducing Cp, increases ship response; the 
greater effect being caused by varying Cp. The JLM designs 
bear this out. Ship JLM-1 is the worst performer, relative 
to the others, since it has the lowest Cp. Ship number 
three is best to 26 knots due to its high Cp. Above 25 
knots the curves of JLM's 3 and 4 cross giving .4 the advan- 
tage at 30 knots. The same thing happens between JLM's 1 
and 5. The crossing is probably attributable to the larger 
displacements of ships 3 and 5. From this one could say 
the ship with the smaller displacement has better seakeep- 
ability at hizsh speeds, but it will not necessarily meét:the 
same military requirements as the’ larger snip. ids 
systematic study of seakeeping trends by varying parameters 
is recommended in the next section. 

The analysis of the 8 feasibility ships were kept 
separate from the 5 JLM ships. The curves of r.ms,. relative 


motion vs. ship speed (Figures VIII through XIII) showed 
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a definite difference in r.m.sS. response between the two 
groups, especially at the lower significant wave heights. 
JLM-2 and feasibility design 7 curves were expected to 
appear fairly close together since their characteristics 
have only small differences. They were not close together 
and had different shapes. This led the author to keep 

the two groups separate pending further study. Inputs to 
the seakeeping subroutine appeared correct for all 13 ships 
but the comparative snape of the response curves showed an 
unaccountable difference. This type of study is beyond the 
scope of this thesis but definitely should be carried out 
before the program is used. 

The performance curves (Figures V and VI) are an 
excellent way to compare a varied number of ship designs. 
The individual differences in performance due to variation 
of parameters and dimensions are weighted and compared 
statistically so that the final result is a direct compar- 
ison of performance only, based on the chosen criterion. 

The oversll result of this thesis then is the establish- 
ment of a logical metnod of evaluating the seakeeping performn- 
ance of a variety of shins during the feasibility design 
Stage from predicted performance curves based upon selected 


criteria, and applicable to computer calculation. 





CONCLUSION 

The main objective of this thesis has been the 
establishment of a method whereby an analyst-designer 
would be able to assess the seakeeping performance of 
many ship designs while still in the feasibility stage 
using currently available computer programs and techniques. 
This has been done. The resulting performance curves can 
be used directly by the designer to evaluate the probable 
performance of a large number of ships. Comparisons can 
be made and decisions arrived at which would ollow trade-offs 
in performance, costs, and meeting military requirements. 

No attempt was made by the author to establish trends 
in seakeepings performance by varying ship characteristics. 
It has been swown 2], (ley leoir that increasing the length 
to draft ratio and beam to draft ratio tend to improve response. 
An increase in waterplane inertia by increasing beam, holding 
Cx and draft constant, leads to V-form section shapes. V-forms 
perform very favorably in a seaway as was demonstrated in [2] 
and snown here by the hizh performance of JLM-4 in Figure VI.* 

Ships of longer length hold a definite advantage in 
reduction of ship motions. However, draft is also beneficial 
wnen evaluating certain criteria, such as sonar dome immersion. 
Deeper draft vessels tend to have less overall bow relative 
motion, Keeping the sonar dome immersed. This fact also would 
help to téduce slamming and water over the deck. Altering a 
particular hull parameter in order to improve a desisn's seak-e 
keepability necessitates a careful review of the effects on 
other parameters which may be changed by the alteration, causing 
adverse effects, possibly offsetting the desired inprovement. 


The technique presented permits this to be done. 
a #See body plan in Confidential REESE 





DISCUSSION AND nis;COMMENDATIONS 


The predicted performance curves arrived at as a 
result of this thesis apply only to the assumed criterion 
of sonar dome immersion. A more detailed seakeeping an- 
alysis of a group of shivos would include performance curves 
based on other criteria. The present program computes 
statistical values for acceleration and velocity of bow 
meeacive vertical motion, stern acceleration and pitch 
anzle. Any criterion, once established, that would make 
use of tne above statistical values could be investigated 
using the proposed metnod. For example, an analysis in- 
volvins slamming would make use of the statistical values 
computed for bow relative vertical acceleration. The pro- 
Meeeo method could be used after establishing a criterion 
Timitine the relative acceleration to a value wnich would 
Moo produce slams. For further discussions of particular 
proposed seakeepinz criteria see(2} andfi8l. 

The predicted performance calculations are not pro- 
sra*med. The complete computer study would involve pro- 
gramming these calculations as a subroutine at which time 
the selected seakeepins criteria may be included. 

The present program is written for only ahead, irrez- 
ular, long crested seas. For complete analysis the calcu- 
lation of snip motions at all headinzs and the ability to 


analyse snip resvonses in short-crested seas should be made 


in 
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as available subroutines to be called as the user desires. 
One must not forget, however, that the application 
of a seakeeping study at the feasibility stage requires 
rapid analysis of a great number of computer designs. 
Any lengthening of the program should be considered care~ 
fully. The program as used herein is already large and 
considered by the author to need careful revision before 
it can be included as a subroutine to a computer feasibility 
study. (To compute the analysis of one ship used in this 
study required a minimum of one hour; 5 minutes on the 
IBM7090 computer, 5 minutes to print the output, and 50 
minutes to tabulate and calculate the predicted performance 
percentages.) Calculations using fewer ship speeds: and 
larger increments of frequency of encounter will reduce 
computer time, but the major reduction would have to be 
program revision based upon future work in simplifying the 
ship motion computations. 
Due to the number of ship characteristics required 
as input to the program, the author believes a study of the 
effects of variation of hull parameters on ship response might 
well be carried out. A variation of any one or a combin- 
ation of the input parameters, Cx, Cp, B, H, La. LCG and 
radius of gyration, carried out in a systematic manner would 
SBRAbL i oh definite trends in overall ship response. The 


naval architect would then be able to focus his attention 
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on the combinations of ship characteristics which would 
result ina ship with good seakeeping qualities. 

The control program written by the author can be 
used with any seakeeping program where values of bean, 
draft and section area coefficients are used to calculate 
the snip responses. It can be detached from the seakeeping 
subroutine used by the author and by changing a few state- 
ments attached to any new or improved seakeeping program. 

It can also be used alone to provide a designer with points 
on a sectional area curve or points on a load water line 
curve. These curves will, however, form a family with the 
JLM designs. 

The complete program, called KSEA, is not in its 
ideal form. Durinsz this study the subroutines would not 
directly accept the data computed by the control progran. 
Cards had to be punched by the control program then fed 
to the subroutines even though the data was identical to 
that stored in the memory of the computer. This should 
be corrected so the program will be continuous. The sub- 
routines should be revised to compute motion and its deriva- 
tives at the forward perpendicular as well as the slamming 
station. The probability analysis of Solution ~- Part IIT, 
should be prozrarmed and made a subroutine of KSEA, possibly 
usinz a plottinz routine to plot performance curves. These 


additions and revisions would make KSEA a functional design 
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computer program accepting as input, Cp, Cx, Bx, Hx,A, 


and L from the feasibility study and giving the designer 
a performance curve, based upon the selected eri ter lence 


he can use directly in his evaluation of the designs. 
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APPENDIX A 


Sea Spectrum Calculations: see also 11 and 18 
The original program assumed a unidirectional sea- 


way represented by the Neumann energy spectrum 


~ KNOB 
(We) 


! ¢ 


me cverms of frequency of encounter, We , Where mm is Wind 
speed, G is wave frequency, and V is ship speed. 
The I.T.T.C spectrum also assumes a unidirectional 
seaway and is represented by 
-33.5@ 
i - 0.00810 Y” exp ee | 
i Ww WG EEE 
mameeerms Of frequency pie encounter, We , where Hy, is 
significant wave height as given in Figure IU. 
The response spectrum is the result of multiplying 
: the sea spectrum by the response amplitude operator. The 
mean-square value of record of the response is the inte- 
grated area under the response spectrum. The root-mean-) 
Square of record is the square root of tne above area. 
The following statistical values can be expressed 
in terms of the root-mean-square value: 


Average Amplitude .886 1 RMS 





Significant Amplitude (average one-third hignest) 1.416 ¥RMS 


Average of one-tenth hiznest amplitudes 1.30 RMS 


% 
¢ 
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These three statistical values are computed for each of 


the five cases: 


1. vertical acceleration at slamming station 

2. relative vertical motion at tne slamming station 
3. relative vertical velocity at the slamming station 
4, relative vertical motion at the propeller station 


5S.  piteh Aanaie 


The slamming station is located .1L aft of F.P. 
and the vropeller station is A.P. 

A typical comparison of the two spectra is given in 
Figure VI. It will be noted that the Neumann Specrum is 
more severe and hence sives larger values of RMS Relative 
Motion. The I.T.T.C. Conference in December 1966 chose 
its spectral formulation based on the work of Pierson and 

Moskowitz, giving a closer approximation to actual seaway 
Comgitions. 

The program, as it presently exists, computed the 
output for this thesis in five minutes per ship case, for 
Bmp input of five snip speeds and fifteen frequencies of 
encounter. The eight feasibility ships studied herein re- 
quired a total of 39 minutes of computer time. 

A more detailed discussion of tne seakeeping program 


fil be four@ in [18]}"na [2]. 














APPENDIX B 
Input and Output Data for Computer Program. 


Input Data. 

Unless stated, any consistant system of units may 
be used in the program, the decision being made by the user 
in the fourth data card. The program is written so that 


computations can be performed on any number of ships. 


Data card 1, Format (4F 10.4) 
Bx - Maximum beam 
Hx » Maximum draft 
Cp - Prismatic coefficient 


Cx = Maximum Section coefficient 


The following data cards are identical to the MIT Seakeeping 


program. (1 8] (2] 


Data card 2, Format (A6, I2, 10A6) 

Identification card. Information on this card will 
be printed at the top of each page. ZK1il must appear in 
columns 5-8, otherwise the user is free to punch any other 


information on the remaining columns. 


Data card 3, Format (5110) 
This card controls the form.of output. There are 


” 
five options. For statistical values to be printed, a1 


Lyk 


she ok 


on a 


oe — eas 











5, =. 


in column 50 instructs the computer to print them. (Output 


Form 5) 


Data card 4, Format (110, 410.4) 


Five values should be given on this card: 


N = the number of stations. This number will always 


be 20 if the entire program is run. 


L - Length between perpendiculars 


Y - specific weight of water - 0.0286 for salt water 


when displacement is in tons. 


& =- gravitational acceleration, 32.2 ft/sec“. 


2 


Av - Displacement of ship in tons. 


he 


Data card 5, Format (2F 20.4) 


This data card should give two quantities: 


1) Radius of gyration (.25L if not known or assumed) 


2) Distance from F.P. of center of gravity, (.5L if 


not known or assumed) 


Data card 6, Format (7F 10.4) 


This data card should give the seven quantities: 


h 


@ min 
GO) max 
Sw 

V min 


VV max 


6V 


wave amplitude (must be unity for statistical 
calculations) 

Smallest frequency of encounter, 0. 

Largest frequency of encounter 

Increment of frequency of encounter 

Smallest ship speed - ft/sec 

Largest ship speed - ft/sec 


Increment of ship speed - ft/sec 


_ - 
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Data card 7, Format (3F 10.4) 

The following information should be punched for 
Neumann Spectrum calculations: 

\Wmin - Smallest wind speed, knots 

Wmax - Largest wind speed, knots 


& W - Increment of wind speed, knots 


The following information should be punched for 
met. t.Cz Spectrum: 
| | min - Minimum significant wave height - ft. 
|| max - Maximum significant wave height - ft. 


SH - Increment of significant wave height - ft. 


mewe: There are oe Te of the above program. One 
uses the Neumann Spectrum, the second utilizes 1.T.TeC. 
spectrum. Both decks are in binary lanzueage (macnine 
lanzuage). 

Tne data cards must be followed by a blank card which 


will end the computer run. 
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MEP UNDAx 6. Sample analysis calculation using proposed method. 


This appendix was added in order to demonstrate the 
entire process of going from hull dimensions to performance 
curves by the use of a computer program. Feasibility design 
number 6 will be used for demonstration purposes. The assump- 


Tions are identical to those stated in Solution, Part Il. 


The KSEA program first reads as input the values of 
Bx = 48.51; Hx = 17.08; Cp = 0.60; Cx = 0.80, taken from 
the feasibility output (Table C-1). These values are used 
in the control program (Table C-2) to compute beam, section 
area coefficient and draft at each station 0-20 (designated 
1-21 for computer application). These values are then stored 
in memory cells designated DRAFT, SECOE and BSTAR respectively. 
They are then called by subroutine SEAKEP and are 
printed, along with other input data listed in Appendix A. 
Table C-III snows the actual computer printout of input data. 
The left column is the waterline curve of beam, the middle 
column is the section area coefficient, the tnird column 
is station draft. These values were computed by the main 
program. Other data is placed in the computer on punch 
cards using the formats listed in Appendix A. The computer 
then calculates statistical responses for each speed selected 
at each wave height over the complete range of frequency of 


encounter. The sample output for a ship speed of 15 knots 
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in a Wave heicht of 34.6 feet is shown in Table C-IV. The 
circled value of r.m.s. relative vertical motiom dn Taba 
C-IV was tabulated at 15 knots over the range of sizgnifi- 
cant wave heights, 14.5 - 48.0 feet, and used in line 3 

of Table V. This was done at each speed for each of the 

i3 ships. The r.m.s. value for significant wave heights 
below 14.5 feet was extrapolated from the r.m.s. vs. Wave 
height curves at constant speed (Figures VIII through XIII). 
The tabular method was then carried out as described in 
Table V, resulting in the figure of 88.1% of time that the 


sonar Will be operable at 15 knots for design number 6. 
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4 Le. 12 17654 2028 ~ 
OM 76 an 34ed4 3253 
Fa3 Pilg 2MeSt.. 27.62: 
KG 18246 L7e95 Te 7.2 
GM elt Ge 9b he 90 
Cv 2024 2002 -" 4.90 
B/H 2.87 2-83 "4 See es" 
VLR 1. 44 Cy ee See. 
FR 2.69 6069 Ge69 
SHP L402 9°7O16 924636 
PI. . alt aga STOLL 9462. 
SHPC 6564. a5979. 158226. 
ELKW 24046 2385. (2354. 
COMPLEMENT st 
COMP @ side 235.6 2aae? 
OFFR aoe shed Gor “19 
CREW 2. & aie. Zire... 
SPACE : — 1% » 
THY SLEGnT. 400714, “S2895-.  449126.- 
SSv GHRES Set fe, 597196! 58656. 
AAH ZLNGZ He CEG * 2uOoecee 206 38. 
AASS 7 he We 6977. C9VESe 6957. 
x AAH Rs oe 8% 0. 
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TABLE C-Tl - List ef Contre! Program, KSEA 


C PROGRAM FUR INTERPOLATING AREA, BEAM,DRAFT AT STATION FOR ENTRANCE 
C INTO A SEAKEEPING PROGRAMe 
C MODELS JLM SERIES 1 


500 
900 


700 


50 | 


Teast ~- EFN SOURCE STATEMENT -— IFn(S) - 


CURVES SOLVED ARE BASED ON SEAHAWK 
THRU Se 

DEMENSTON Y(81).TR(6)-TI(6),ADDA(4) -BEEB( 4) -CGGC(4).UR(6) -UI(6)> 
LOMASS(21) eQUANT (21) oe SKLAM(21) eBSTARC21 ) eCXFST( 21) e¢SKXFST(21) » 
2CTFST(C 21) eSTFST(21) oe XTC 289 -DI XC 21) -ENOXI(21),D0RAFT( 21) OWEIGH( 21 Jeo 
3SECOE(21),ABAR(21),SPECM(B1 ) s VOMEGA( 81) »sRAOSDODO(61) »-RAORB( Bl). 
S4RADROB(81) eRAORP( B81) RAOTHE (81) -RSDSD00(81)-RSDRA( B61) -RSORDB(SI),. 





SRSORP(S1) -RSDTHE( SL) . TITLE 12) 
6 .SP(18), OM( 40). ZN(640)e DE(40)e IN(40). EP(40)+ OME(4) 
7 sAN( 21) sBN(6 21) HNO2E) »CAS(21 De AS( 21) eBS(24) HS(21) 

COMMON Y Pee » SYMPS  . DXI » ADOA BEEB 
COMMON CGSC >» ZREAL »« ZIMAG » FREAL »« TIMAG ZNULL 
COMMON TNULL » DELTA ., EPSIL » TR - ft v 
COMMON ODMASS , QUANT ¢« OMEGAE « SKLAM « KRIT RO 
CoMMON GRAV e BSTAR eo CREST e SXKFST e ALPHA SHNULL 
COMMON xX I » BETA eo BMNULL »« GAMMA » OX M 
COMMON WA ° WAVEN e¢ CW © ENOXS +» SIGMA TAU 
COMMON FNULL . EMNULL » DRAFT » OWEIGH , SECOE TMASS 
COMMON WN » UR . ul e ABAR » PY SHREAL 
COMMON SHIMAG « BMREAL » BMIMAG »« YNERT « VOITH RIZA 
COMMON OMEGA ¢ MM » NN » H . KK Ww 
COMMON VOMEGA , SPECM »« RAOSDD « RAGRB .« RAORDS RAURP 
COMMON RAOTHE . RSDRB «+ RSDSOD ». RSDRUDB . RSDRF RSDTHE 
COMMON RSTSOD ¢« RSTRP e RSTROD , RSTRE eo RSTTHE HELP! 
COMMON HELP2 .»« HELP3 »« HELPS » HELPS . AVESOD AVERB 
COMMON AVERDB . AVERP , AVETHE . SIGSOD , SIGRB SIGRDB 
COMMON SIGRP , SIGTHE »« AVHSDD « AVHRB » AVHRDH AVHRP 
COMMON AVHTHE . CONST » GSQUAR . OFENFR 

FQUI VALENCE (DRAFT.HS) .(BSTAR. BS)» ( SECOE.CAS) 

READ(S,900) Bk ehHX eCPeCX 

FORMAT(4F 1004) 

IF(BX) 6002600,700 

AX=CX#BXeHX 

CP2=CP¥e¥*2 ~ 

BN(1L) = Oc387108— 14256049eCP# 120274998CP2 

BN(2) =—-Se3Bl17994+17. 771 660*CP-142416660¥CP2 

BN(3) =-040416799+ 1eV13IIZECPH 0. 166066*CP2 

BNL 4) =+003941994 t.1449998CP4 0.250000#CP2 

BN(S) =-0.3175994 1.t 7E6GE6*CPH 0.3333334CP2 

BN(G) = 461035199-14.111660*CP4! 3.916660¥"CP2 

BN( 7) =-0.0641994 164366664%CP-0. 1666664%CP2 

BNC8) =-06612999+ 36859999#CP— 2.499999*CP2 

BN(9) =-0.7331994 S4.819998CP- 3464999998CP2 

ON (10)=-00e0837994¢ 3.1133334CP=— 24333333*CP2 

UNC LL) 0c6376004 1LeO0849498CP— 0.6749999"%CP2 a 
ENC 12)=02922600+ 043233334CP-0.333333*CP2 

BANC 13)= 069072004 046304999*CP—- 02249999*%CP2 

ENC 14)= 069321994 0.2550004CP—0.250000#CP2 

BNC I1S)= 025440004 146233333*CP= 0.833333*CP2 

BN(16)= 0.6.372599+¢ 1e415000%CP— 0.750000*CP2 b 
BNC 17)= 02565000+ 0.30683338CP+ 0641 6666¥CP2 i 
BNCLB)= 22256599- 5.951666*CP+ 52916666*CP2 

BNC IGD)= 02232599 0+.5400008CP+ 0.499999*CP2 

BN(20)= 0.0805994+ 0.831666*CP+ 0.08333*CP2 





SO 


300 
903 


400 


600 


Test - 


BN(2l)=1-502599- 


AN( 1) =-O0-1091 22+ 
AN(2) = 02040989- 
AN(3) = 0-0862416- 
AN(4) 00125259- 
AN(S) = 0.2272381- 
AN(6) = 0.2£542730- 
AN(7) = 02440767— 
AN(8) =-02482566+ 
AN(9) =-0.273566+ 
AN(IO)= O2131915+ 
ANC LT)= 126052243- 
AN(12)= 02e95074%6+ 


AN(13)=-0-022476+4 
AN(14)=-12056005+4 
ANC 15 )=-2-e904320+ 
AN( 16)=-12106661+4 
AN(C17)=-0.073434+ 
ANC 138)= 02158555- 
AN(C19)= 02607695- 
AN(20)=-0-039063-— 
AN(21)=-0408281 5+ 
DQ 5O IT=1e15 
HNC [T )=LeO 


m. 


EFN 


02-5164 02*"CP=- 
0e2555208CP4F 
0-2-4534 39#*CPe 
0-540 741 *CP+ 
0.8571 43ecP+ 
1-498412*CP+t 
0.-6947098#CP+ 
3e021957*¢CP- 
2e821957#*CP- 
2e067989*CP—- 
0-632804*CP+ 
Oel 76984€CPH= 
3e1214284*CP- 
6-147089*CP—- 


SOURCE STATEMENT 


3255999€CP4# 2.2.999999*CP2 


0+ 4232808#CP2 
0.6867831*CP2 
1.864021*CP2 
164814381 *CP2 
16904761 *CP2 
20539628¥*CP2 
16 7989418CP2 
1-534391 @€CP2 
10534391 *CP2 
1e216931#CP2 
0.846561 *CP2 
0.158730€CP2 
2. 380952*CP2 
4 0656084*CP2 


13 «Si 9S704CP-8.677248*CHP2] 


4.286243*CP- 
Oel7BO042S*CHP + 
O.-8777778CPF 
229021 16*CP+ 
Oe LSISS87*CPe 
024351 85%CP-— 


22010582*CP2 
LeS534391*CP2 
2e222222¢CP2 
3$e280423*CP2 
Oe 793651 4*CP2 
02 370370¥CP2Q 


IF (CP eL Te eSE HN 15)=02-0434+1.650*CP 

HNO 16) =~-1LeO72399 C4 -731LEE66*CP—2-416666#CP2 
HN( 1 7)==+026337994+246966664CP— 0e 6666664CP2 
HN (18) =-— 1216299944, 062499#CH-1-28749998CP2 
HN (£9 )=—-OeM1979G F121 7499"CP 4+ O00 3750008CP2 
HN ( 20 ) =~ 1460439945635 6666#CP=- 3.666666FCP2] 
HN(21)=02463399-1.160833#CP4+0-9583334CP2 


DC300 IT=1e21 

AS( IT )=AN( IT) AX 
BS( TI )=BNC I) *BX 
HS( 1 )=HNC ET) #HX 


CAS(CI)=ASCTIIS(IBSC IT) FHS(T)) 


READ (5-903) 
FORMAT (F100) 


CAS (1) 


PRINT 400e6CASCIE) BSC EVPHSCI)-CAS(CI),T=1221) 
FORMAT(4(F100e4.5X)) 


CALL SEAKEP 
GO TO 500 
oP 

END 
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Calculation of percentage of time that bow sonar would be 
operable ~- Feasibility Design 6 at 15 knots in assumed 
Wave distribution. 














3 TeMeSe 


(a) 


(b) 


ec ) 


(ad) 


eee? 14.5 14.5 Zine 279 34.6 414.3 48.0 
r 


2 Wave frequency 0.737 0.156 0.069 0.031 0.0045 0.0020 0.0005 


+ Product 1.695 0.839 0.559 0.320 0.056 OvGad 0.009 
lines 2 x 3 
Note (b)-—~ = .16(3.509) = .561 Total line 4 = 3.509 


p xe OeiBn-dine 3 6.918 1.426 -3.438 -7.403 -11.304 -16.206 -21.68 


5 Probability of 1.0000 .9230 O0-5000G 070000 “O7e000 0.0000 0.000C 
not exceeding 
Tetlle O= eei8 


milines 2 x 6 0.737 0.144 0.0 0.0 513 0.0 0.0 


TABLE C-V 


value (2.30) 5.380 8.108 10.332 12.521 15.269 ieee 
relative motion 





Total Line 7 = .881 


It is assumed that at any particular time the bow sonar 
is considered operable if r.m.s. < 6.18 ft. This value 
is based on a limiting value of (H = 1) or 16.08 ft where 
H=17.08 ft. (From Figure III) 


Standard deviation,o’, assumed to be .16 x total of line 
four or .16 x weizhted mean r.n.s. 


The quantity X , line 5 is tabulated in statistical hand- 
books from which line 6 is obtained. 


The total of line 7 gives the probability that the sonar 


will be operable ina unit time. For this illustration 
88.1 %. This point is then plotted. 
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